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A  spacially  extended  system  of  two-level  a'oms  interacting  v  ::h  a  rcso"?r.t  c'cctror.'.agiicuc  field  exntbiis  cunditio.is  for  wliich 
two  phases,  of  ntgh  and  low  excitation,  coexist  spacially,  in  the  direction  of  propagation  in  the  material.  The  spacial  first-order 
phase  transition  and  its  properties  are  related  to  an  inversion-dependent  renormalization  of  the  resonance  frequency,  which  becomes 
significant  for  a  collection  of  two-level  atoms  with  high  density  and  large  oscillator  strengths. 


The  phenomenon  of  optical  bistability  <OB)  is 
usually  based  upon  two  components:  (a)  The  illu¬ 
mination  of  an  optically  nonlinear  material:  (b)  A 
feedback  mechanism  which  is  provided  by  the  mir¬ 
rors  of  a  htgh-0  cavity  [  1  ].  It  was  first  shown  by 
Bowden  and  Sung  [2]  that  OB  without  external 
feedback  may  occur  for  a  system  composed  of  a  col¬ 
lection  of  two-level  atoms  interacting  via  the  electro¬ 
magnetic  field  and  driven  by  an  externally  applied 
coherent  field.  By  using  a  fully  quantized  many-body 
model  in  the  mean-field  limit,  they  have  shown  [2] 
that  the  role  of  the  feedback  may  be  provided 
intrinsically  by  the  induced  coherent  dipole-dipole 
interactions  among  the  atoms  of  the  material. 

For  sufficiently  low  values  of  the  incident  field,  the 
reaction  field  produced  by  the  induced  dipoles 
opposes  the  externally  applied  field  and  thus  reduces 
the  net  field.  For  sufficiently  large  values  of  the  inci¬ 
dent  field,  the  individual  atom’s  interaction  with  the 
electromagnetic  field  is  influenced  much  more 
strongly  by  the  externally  applied  field  than  the  reac¬ 
tion  field  due  to  the  interatomic  interaction.  The 
incident  field  overwhelms  the  interatomic  interac- 

1  A  National  Research  Council  Research  Associate  Permanent 
home  institution:  Technion  -  Israel  Institute  of  Technology, 
Department  of  Physics.  Technion  City,  32  000,  Haifa,  Israel. 

;  A  National  Research  Council  Research  Associate.  Present 
address:  Department  of  Physics,  Southern  Methodist  Univer¬ 
sity,  Dallas,  TX  75275.  USA 


lion  via  the  electromagnetic  field,  and  thus  the  con¬ 
tribution  of  the  reaction  field  due  to  the  coherent 
dipole-dipole  interactions  is  largely  negated.  This 
occurs  as  a  first-order  phase  transition  far  from  ther¬ 
modynamic  equilibrium. 

Recent  interest  in  intrinsic  optical  bistability  ( IOB ) 
has  stimulated  a  surge  in  theoretical  and  experimen¬ 
tal  work  [3-13].  Particularly  relevant  to  the  work 
presented  h"re  is  the  recent  theoretical  discussion  of 
the  effect  of  the  insertion  of  the  local-field  correction 
in  the  usual  Maxwell-Bioch  formulation  in  causing 
IOB  [12],  and  the  origin  of  the  need  for  that  correc¬ 
tion  in  the  passage  from  microscopic  to  macroscopic 
electrodynamics  [13]  in  semiclassical  theory.  The 
origin  of  IOB  was  identified  in  that  treatment  as  due 
to  an  inversion-dependent  renormalization  of  the 
atomic  transition  frequency  caused  by  coherent 
dipole-dipole  interaction  among  two-level  atoms. 
Previous  quantum  mechanical  treatments  of  the 
problem  [2]  used  a  mean-field  approximation  in  a 
volume  much  smaller  than  a  cubic  wavelength  and 
suffered  from  the  defect  that  propagation  and  retar¬ 
dation  effects  were  neglected.  Although  the  treat¬ 
ments  using  semiclassical  theory  include  these  effects 
within  the  framework  of  the  semiclassical  approxi¬ 
mation.  the  recent  work  was  been  addressed  to  either 
a  propagation  length  much  smaller  than  a  resonance 
wavelength  [12]  or  to  a  volume  containing  a  speci¬ 
fied  ( small )  number  of  atoms  [13]. 
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In  more  recent  works  [14-16]  we  consider  the 
problem  from  a  many-body  quantum  mechanical 
approach  in  which  we  treated  a  large  number  of  dri¬ 
ven,  interacting,  spacially  distributed  two-level 
atoms.  The  intrinsic  many-body  statistics  justify  fac¬ 
torization  in  the  multiplication  of  the  dipole  opera¬ 
tors  [15],  while  propagation  and  retardation  effects 
are  included  and  found  to  be  important  [14,15],  It 
has  been  shown  [  1 4, 1 5  ]  that  interatomic  correlation 
is  built  up  in  the  direction  of  propagation  of  the  inci¬ 
dent  field  and.  that  under  appropriate  conditions  and 
limits,  the  usual  Maxwell-Bloch  equations  are  repro¬ 
duced.  However,  for  high  densities  and  large  oscilla¬ 
tor  strengths,  frequency  renormalization  which  is 
inversion-dependent  appears  as  a  modification  to  the 
usual  Bloch  equation  for  the  complex  atomic  dipole 
operator  [14-16], 

In  this  letter  we  show  that,  due  to  the  nonlinear 
frequency-renormalization  correction,  a  new  kind  of 
spacially  dependent  first-order  phase  transition  may 
occur  in  the  interaction  between  a  coherently  driven, 
interacting  collection  of  spacially  distributed  two- 
level  atoms  of  large  density  and  high  oscillator 
strengths.  We  show  how  the  two  phases  of  high  and 
low  excitation  may  coexist  spacially  in  the  propaga¬ 
tion  direction  in  the  material.  Due.  in  part,  to  this 
longitudinal  spacial  dependence,  the  first-order  phase 
transition  predicted  here  for  IOB  is  fundamentally 
different  from  the  well-known  effect  of  bistability  in 
resonators. 

The  equations  of  motion  derived  by  us  are  given 
as  [14.15] 


d<i7.  >/d/=  —  ;•(  [  < <r.  >  -H  ] 

+  (/i/fi )[E*< <7  .„>  +£<(7  „>  ]  . 

(1) 

d<a  . >/d i-  -  (nFJ2fi){  a .  > 

+  i[J  — «<<r, >]<(;.„ >-yT<ff  .„>  . 

(21 

(\/c)dE/dz=  -  (4rtw/i/c)<cri|,>  . 

(3) 

Here.  E  is  the  slowly-varying  envelope  of  the  electric- 
field  amplitude,  <<7.„>  is  the  expectation  value  of 
the  slowly-varying  complex  atomic  polarization  per 
unit  volume,  and  <<;.->  is  the  expectation  value  of 
the  atomic  inversion  per  unit  volume.  The  equations 
for  Zf*  and  <  er  0  >  are  the  complex  conjugates  of  eqs. 


(3)  and  (2).  respectively.  The  parameter  fi  is  the 
modulus  of  the  matrix  element  of  the  transition 
dipole  moment  of  an  atom,  co  is  the  frequency  of  the 
incident  field  and  n  is  the  number  of  atoms  per  unit 
volume.  The  rates  ;t  and  y  T  are  the  inverses  of  the 
relaxation  times  7j  and  7\  for  the  inversion  of  pop¬ 
ulation  and  the  dipole  moment,  respectively,  and 
J  =  u>-W(i  is  the  deviation  of  the  applied  field  fre¬ 
quency  oi  from  the  atomic  resonance  frequency  can. 
Finally,  t  is  the  frequency  renormalization  constant 
derived  in  our  previous  work  [14.15] 

t  =  lnnpc'/4w'i  .  (4) 

where  fi  is  the  spontaneous  decay  constant 

(l  =  4\fi\2(Ou/ificx  .  (5) 

The  frequency  renormalization  stems  from  coherent 
dipole-dipole  interactions  between  atoms  that  are 
within  a  volume  of  a  cubic  wavelength  [2.12-1 5). 

In  the  fullv  quantum  mechanical  model,  y,  =/f. 
y  ,=/?/2.  but  by  considering  additional  homogene¬ 
ous  broadening,  we  may  treat  ■/,  and  as  empirical 
constants.  The  usual  Maxwell-Bloch  equations  [17] 
are  obtained  from  eqs.  ( 1 )-( 3 )  in  the  limit  e  ->0.  The 
frequency  renormalization  correction  proportional  to 
(  becomes  important  only  for  a  two-level  system  with 
high  oscillator  strengths  and  a  large  density  .  Such  a 
system  can  be  realized  in  experiments  with  Rydberg 
atoms.  While  the  analysis  presented  here  is  suitable 
mainly  for  an  atomic  system,  it  suggests  that  similar 
effects  may  be  obtained  by  using  crystals  with  a  high 
density  of  bound  excitons  having  large  oscillator 
strengths,  which  may  be  realized  in  CdS  for  bound  I: 
excitons  [18], 

Under  steady-state  conditions  we  get.  from  eqs. 


(1 ) — <  3 ) ; 

<(T.„>  =nE{a:  }/2fi[iJ-  y,  -  ie  <  a.  }  ]  .  (6) 

7i.[<0  +  ■  ]{7f  +  [J-tztr.)  ]-  ’ 

—  -/i  /<<?.'/  .  ( 7  ) 

<5//df=  -x,  [  I  +  <<t->]  ,  (8) 

where  /  is  the  normalized  intensity  of  radiation. 
!=\nE/h\\  (9) 

and 
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Fig.  1.  The  atomic  inversion  of  population  <<t.>  is  calculated  as 
a  function  of  the  normalized  distance  2  [defined  in  units  of  s  1 
by  eq.  (10)]  using  the  normalized  parameters:  yL  =  2,  y,=  l. 
d=  -  4  ( in  units  of  s  1 ).  The  intensity  of  the  incident  radiation 
/  [defined  in  units  of  s  ;  by  eq.  (9)]  is  given  by  /(r  =  0)  =  100. 
Curves  ( a )  and  ( c )  are  calculated  according  to  eqs.  ( 7 )  and  ( 8 ) 
with  the  parameter  t  =  20  (s  1 )  using,  respectively,  the  upper  and 
lower  branch  solution  of  the  cubic  equation.  Curve  (b)  is  calcu¬ 
lated  according  to  eqs.  (11)  and  (12)  with  parameter  e  =  0  (the 
regular  Maxwell-Bloch  equations). 

z=(4nn2(on/fic)z .  (10) 

The  normalized  distance  z  and  the  constants  >'T,  yL, 
J  and  e  are  given  in  units  of  s  1  while  /  is  given  in 
units  of  s 

Eq.  ( 7)  is  a  cubic  equation  leading  to  intrinsic  local 
bistability  because  of  the  nonlinear  term  which  is 
proportional  to  «.  Eqs.  (7)  and  (8)  can  be  numeri¬ 
cally  integrated  to  give  the  spacial  distribution  of 
excitation  in  the  extended  medium,  given  /(z  =  0). 

For  the  usual  Maxwell-Bloch  equations  ( t  =  0 ),  we 
get  simple  analytical  solutions: 

<.o-.>  =  -VL(Y\+f)Hyi.(y\+A2)  +  yTl]  (11) 

representing  the  usual  saturable  nonlinearity,  and 

*  Myf+d2)ln[/(z)//(0)] 

+  yT[/(z)-/(0)]  =  -yLyTz,  (12) 

which  yields  only  one  solution  (no  bistability)  for  I 
t  at  each  z  [  1  ] . 

In  fig.  1  are  shown  results  of  some  of  our  calcula¬ 
tions  for  <  <T->  as  a  function  of  z.  We  have  taken  yL  =  2 
so  that  all  rates,  including  z,  are  in  units  of  (yL/ 2), 
while  /  is  in  units  of  (yJ2)2.  Ml  the  curve*  in  fig  1 


of  <cr_->  as  a  function  of  z  have  been  calculated  by 
using  the  parameters  yL=2,  yT=l,  J=— 4. 
/( z— 0 )  =  1 00.  Curves  ( a )  and  ( c )  have  been  calcu¬ 
lated  according  to  eqs.  ( 7 )  and  ( 8 )  with  the  param¬ 
eter  e  =  20,  using,  respectively,  the  upper  and  lower 
branch  solutions  of  the  cubic  equation  [eq.  (7)].  For 
comparison  with  the  ordinary'  Maxwell-Bloch  equa¬ 
tions,  we  describe  in  curve  ( b )  the  calculation  of  /  as 
a  function  of  z according  to  eqs.  (11)  and  (12).  ( i.e.. 
for  t  =  0 ) . 

By  starting  with  a  high  intensity  of  the  incident 
electromagnetic  field  and  reducing  its  strength,  we 
expect  <<T->  as  a  function  of  z  to  follow  the  upper- 
branch  solution,  of  the  bistable  region,  as  demon¬ 
strated  in  fig.  1  by  curve  (a).  By  starting  with  a  low 
intensity  for  the  incident  field  and  increasing  its 
strength,  we  expect  <a.)  as  a  function  of  z.  for  the 
same  set  of  parameters,  to  follow  the  lower  branch 
solution,  as  demonstrated  in  fig.  1  by  curve  (c).  A 
spacial  first-order  phase  transition  will  occur  for  the 
two  cases  at  different  distances  along  the  long  sam¬ 
ple,  for  the  same  incident  intensity.  In  order  that  the 
effect  will  be  pronounced,  the  value  of  e  should  be 
large  (high  density  and  large  oscillator  strengths). 
Also,  in  order  to  get  the  phase  transition,  the  fre¬ 
quency  of  the  coherent  radiation  should  be  lower  than 
the  atomic  transition  frequency  and  of  such  a  value 
that  the  frequency  renormalization  will  bring  it  in  and 
out  of  resonance  as  a  function  of  the  value  of  <a.) 
(as  demonstrated  in  fig.  1  by  using  the  values  J  =  -4. 
€  =  20). 

To  obtain  a  good  condition  for  intrinsic  bistabil¬ 
ity,  the  parameter  t  should  fulfill  the  condition 
e/yT>  6.  where  yT  is  the  relaxation  rate  for  the  com¬ 
plex  dipole.  Assuming  for  the  pure  quantum  case 
yT=/?/2  where  /?  is  the  spontaneous  decay  rate,  we 
have  f//?>  3  and  we  find  that  we  need  approximately 
a  few  hundred  (or  more)  atoms  within  a  volume  of 
a  cubic  wavelength  in  order  that  the  effects  will  be 
significant.  If  we  introduce  additional  homogeneous 
broadening,  which  increases  the  value  of  yT,  the  den¬ 
sity  should  be  increased  by  an  additional  factor  2yr/fi. 
The  present  effects  are  therefore  important  only  for 
high  densities  and  large  oscillator  strengths.  This  is 
the  main  reason  why  in  ordinary  transmissivity 
experiments,  these  effects  are  negligible  and  the  usual 
use  of  the  Maxwell-Bloch  equations  ( e  =  0 )  is  justi¬ 
fied  At«e.  C”“  sb^'ild  note  that  the  longitudinal  spa- 
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cial  first -order  phase  transition  is  predicted  here  as 
an  internal  effect  in  the  material  for  the  inversion  of 
population  along  the  sample,  for  certain  critical 
parameters.  The  observation  of  this  effect  needs, 
therefore,  the  use  of  transverse  probe  fields  by  which 
local  inversions  of  populations  are  detected. 

In  the  present  analysis  we  have  neglected  the  effect 
of  output  boundary  conditions,  which  is  a  quite  good 
approximation  for  a  very  long  sample.  Also,  we  have 
assumed  in  our  model  a  homogeneous  broadening. 
The  spacial  first-order  phase  transition  obtained  in 
the  present  theoretical  analysis  is  a  fundamentally 
new  effect  in  quantum  optics.  We  predict  here  the 
possibility  for  its  observation  in  experiments,  where 
the  inversion  of  population  <a_.>  can  be  measured 
as  a  function  of  z  by  using  a  probe  electromagnetic 
field  propagating  in  the  transverse  direction  of  a  long 
sample. 
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